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We show that the observed rotation curves of spiral galaxies constrain the sound speed of the dark 
matter to be c s < 10 _4 c, where c is the speed of light in vacuum. Using the Modified Chaplygin Gas 
as a representative example of a class of unified dark energy models incorporating an effective dark 
matter component with a non-zero sound speed, we determine the most stringent constraint to date 
on the value of the constant contribution to the equation of state parameter in this class of models. 

Finally, we explain the reason why previous constraints using the Cosmic Microwave Background 
and Baryonic Acoustic Oscillations were not as competitive as the one presented in this paper and 
discuss the limitations of the recently proposed Extended Chaplygin Gas. 


I. INTRODUCTION 

Increasingly precise cosmological observations M re¬ 
main consistent with the standard ACDM model. In this 
model it is assumed that general relativity provides an ac¬ 
curate description of gravity on cosmological scales, and 
that the two main constituents of the Universe are a cos¬ 
mological constant A and a non-relativistic Dark Matter 
(DM) component. While a cosmological constant, or a 
more general Dark Energy (DE) form, is necessary in 
the context of general relativity to account for the cur¬ 
rent acceleration of the Universe, the DM component is 
required in order to explain the observed dynamics of 
cosmological perturbations over a wide range of scales. 
Despite the simplicity of the cosmological constant, there 
is presently no satisfactory explanation for its tiny energy 
density, which favours dynamical DE or modified gravity 
as a better motivated explanation for the acceleration 
of the Universe (see, for example [Ml and references 
therein). In fact there is no alternative to dynamical DE 
or modified gravity if the vacuum energy is screened, and 
prevented from acting as a gravitational source (see, e.g., 

BE!). 

Another (albeit equivalent) interpretation of the 
ACDM model relies on a single perfect fluid with con¬ 
stant negative pressure. Hence, ACDM may be regarded 
as the simplest Unified Dark Energy (UDE) model, in 
which DM and DE are taken as different manifestations 
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of a single unified DE fluid [H|. A broad class of UDE 
models, known as the Generalized Chaplygin Gas (GCG) 
0, includes the original Chaplygin Gas 0, 0 and 
ACDM as particular models. The GCG has been claimed 
to be essentially ruled out due to the late time oscillations 
or the exponential blow up of the DM power spectrum 
predicted using linear perturbation theory [0 (except 
for a tiny region of parameter space very close to the 
ACDM limit). However, it has recently been show that 
the GCG may be consistent with current observational 
constraints, over a wide region of parameter space, as¬ 
suming that there is a sufficiently high level of non-linear 
clustering on small scales 0 (see also 0-0). 

The Modified Chaplygin Gas (MCG) has been pro¬ 
posed as a further generalization of the GCG model 24], 
This model incorporates an extra parameter which is as¬ 
sociated to a lower bound on the sound speed of the 
effective DM component. The MCG has been obser- 
vationally constrained using the redshift dependence of 
the apparent magnitude of type la supernovae, the Cos¬ 
mic Microwave Background (CMB) anisotropies, and the 
matter power spectrum, including Baryonic Acoustic Os¬ 
cillations (BAO) |2a-[3J|. Further extensions incorporat¬ 
ing additional terms accounting for a possible growth of 
the comoving Jeans length with cosmic density have also 
been proposed in the literature (see, e.g., 0). 

In this paper we shall constrain the value of the DM 
sound speed using the observed rotation curves of spiral 
galaxies, and use that result to impose stringent limits on 
the MCG class of models. A comparison with previous 
constraints obtained using other observational data and 
a discussion of the limitations of the recently proposed 
Extended Chaplygin Gas (ECG) will also be presented. 
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Throughout this paper we use units in which the speed 
of light in vacuum is c = 1. 

II. MODIFIED CHAPLYGIN GAS 

The MCG can be described as a perfect fluid with an 
energy-momentum tensor 

T^ = (p + p)uV+pg^, (1) 

having a barotropic equation of state parameter given by 


in the flat part of the spiral galaxy rotation curves. Here 
r is the distance from the galactic center, G is the gravi¬ 
tational constant, M(r ) is the total mass inside a sphere 
of radius r centered in r = 0, p{r) is the average den¬ 
sity inside that sphere, and Vf is the value of the circular 
velocity at large distances from the galactic center. 

If the DM has a non-zero sound speed, then its gravi¬ 
tational collapse can only take place on scales larger than 
the Jeans length [Hj 


XJ[DM ] — C S [DM] 



( 6 ) 


w = B — 


A 

Y+ 1 ’ 


( 2 ) 


where A, B and a are constant parameters, w = p/p, p is 
the proper energy density, p is the hydrostatic pressure, 
are the components of the 4-velocity and g are the 
components of the metric tensor. 

For B = 0 the MCG reduces to the GCG, which is 
arguably the simplest non-trivial class of models where 
the role of DM and DE is played by a single dark fluid. 
On the other hand, for B = 0 and a = 0 the MCG is 
completely equivalent to the ACDM model. In this paper 
we shall mainly be interested in the case where B ^ 0, 
in which case the sound speed is given by 
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2 
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dp 
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Since imaginary sound speeds are associated with insta¬ 
bilities, in particular on very small scales, in this paper 
we shall consider model parameters in the range A > 0, 
B > 0 and 0 < a < 1, for which c 2 > 0. Note that 
the value of c 2 is always greater than or equal to B but 
it tends to B for large enough energy densities. In this 
limit the MCG behaves effectively as DM. 


III. OBSERVATIONAL CONSTRAINTS 
A. Galactic rotation curves 


On smaller scales the DM pressure balances gravity giv¬ 
ing rise to stable oscillations. 

Substituting Eq. © into Eq. © one obtains 

Xj[DM] 
r 

DM may only collapse on length scales smaller than r if 
Xj < r or, equivalently, if 


4 c. 


s[DM] 

Vf 
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C S [DM] < 
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In the context of the MCG UDE model, this leads to the 
following constraint 


B < c 


2 

s[DM] 


3 v j 

4 Y ■ 


( 9 ) 


Here, c s [dm] is to be interpreted as the sound speed 
of the MCG in the galactic halo (again note that the 
MCG sound speed squared approaches B at high en¬ 
ergy densities). Taking into account that the typical val¬ 
ues of Vf observed for spiral galaxies are in the range 
100 kms -1 < Vf < 300 kms -1 |36l - l38j . with some galax¬ 
ies having lower maximum velocities, one finally obtains 
the following conservative limit on the value of B 


B < 10" 8 . 


( 10 ) 


The observed rotation curves of spiral galaxies pro¬ 
vide striking evidence for the presence of DM on galactic 
scales (see, e.g., Yi - IYI ). These curves typically flatten 
at large distances from the galactic center, around and 
well beyond the edge of the visible disks, indicating the 
presence of a (nearly) spherically symmetric DM halo. 
In the Newtonian limit, expected to be valid on galactic 
scales, the circular velocity v is given by 


v 2 (r ) 


GM(r) 
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47T 
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( 4 ) 


so that 


Gp{r) 


4tT r 2 


( 5 ) 


B. Other observational constraints 


The constraint given in Eq. m is much more strin¬ 
gent than those obtained using the CMB or BAO (see, 
e.g., [HJ). Here we shall explain why. 

The adiabatic sound speed of the baryon-photon (bj) 
plasma before last scattering is given by (40] 
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where 2 = 1 — 1/a is the redshift, a is the cosmologi¬ 
cal scale factor, = Pi/p c , an ‘i’ represents a partic¬ 
ular energy component, in this case either baryons (6) 
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or photons ( 7 ), p c = 3H 2 /8 ttG is the critical density, 
H = lOOhkms -1 Mpc ” 1 is the Hubble parameter and a 
‘0’ denotes the present time. Using the constraints on the 
values of the cosmological parameters obtained by the 
Planck 2015 team, namely flboh 2 = 0.02230 ± 0.00014, 
Zrec = 1089.90 ±0.23 and h = 0.6774 ±0.0046 0, as well 
as f2 7 o h 2 = 2.47x 10 ” 5 one may estimate the sound speed 
of the baryon photon plasma just before recombination 


C s[b 7 ](~ r 


, 1 / 3 

c - } "3 ( 1 + 4H 




70 
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( 12 ) 


Here z re c denotes the redshift of recombination and a mi¬ 
nus represents an instant immediately before that. The 
position of the first acoustic peak at £ ~ 200 is extremely 
sensitive to the value of c s r 6 7 ](z T . e c-)- 

A non-zero DM sound speed would be responsible for 
additional signatures in the CMB primary anisotropies. 
In order for them to only affect multipoles with £ > 2500 
(£ ~ 2500 is roughly the upper limit of the multipole 
range probed by Planck, and twice that of WMAP), the 
value of the DM sound speed would need to be 

C s[DM]( Z rec -) < ^2500^ C s[b'y]( Z rec-) ~ 1 X 10 , (13) 

or, equivalently, B < 10 ” 3 if the GCG is considered as 
a UDE model. Not too surprisingly, this limit is of the 
same order as that otained in [32j using the WMAP seven 
year CMB data (together with BAO and type la Super¬ 
novae data). Note that the Integrated Sachs-Wolfe effect 
has also a si gnif icant contribution to the CMB constraints 
reported in [32| (see also 0). 

Given that the Planck results constrain the mat¬ 
ter power spectrum on comoving wavenumbers k < 
0.2 Mpc -1 (42|, the current CMB data is only able to con¬ 
strain the DM sound speed at z rec - on comoving scales 
larger than A > 30 Mpc. Since the evolution of the DM 
cosmological comoving Jeans length during the matter 
dominated era is given by 


A' 


J[DM} 


^ J[DM] 


Cs[DM] \ —1/2 

°c 1/2ffl OC C s[DM] (l + z) ' 


(14) 


it is possible to improve the CMB constraints on the value 
of a constant cA DM , (or B) by a factor of approximately 
1 ± z rec ~ 10 3 using the observed matter power spectrum 
at z = 0 for comoving wavenumbers k < 0.2 Mpc -1 . This 
is consistent with the constraint on the GCG class of 
models reported in [2Sj (B < 10” 6 ). On the other hand, 
Lyman-alpha constraints on the matter power spectrum 
on comoving wavenumbers up to k ~ 2 Mpc ” 1 10 
are expected to lead to even tighter limits on the value of 
B } only slightly less constraining than the ones obtained 
in the present paper using the galactic rotation curves. 


High redshift constraints on the value of the DM sound 
speed may be more effective than low redshift ones in 
models where the Jeans length grows with redshift. An 
example of a class of models where this can be realized 
is the ECG, in which the equation of state parameter is 
given by |35| 


» = AM"?”- 1 --A - as) 

n ' 

Here, n > 0 are integers and B n are real constants. For 
the sake of simplicity let us assume that A = 0, B 2 > 0 
and B n = 0 for n ^ 2. In this case, the ECG sound speed 
is given by 

c 2 s = 2 B 2 p oc (1 ± z) 3 , (16) 


implying that the cosmological comoving Jeans length is 

A / 2 


A 1 


J[DM] 


^ J[DM] 
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during the matter dominated era. This leads to a value 
of A j[ DM ] ( Zrec -) which is larger by a factor of more than 
10' than the value at z = 0. As a consequence, for n > 2 
CMB constraints on the ECG class of models effectively 
rule nearly all the available parameter space, except for 
a very small region with B n ~ 0. 


IV. CONCLUSIONS 


In this paper we have determined the most stringent 
constraint to date on the effective DM sound speed using 
the observed circular velocity of spiral galaxies at large 
distances from the galactic centre. Our results were then 
used to show that B < 10 -8 if the MCG is to be regarded 
as a UDE candidate. We have compared this constraint 
with those obtained by various authors using other ob¬ 
servational data, explaining the reason for the significant 
improvement obtained in this paper. Finally, we have dis¬ 
cussed the case of the ECG, as an example of a class of 
models where the Jeans length may increase significantly 
with redshift, showing that, in this case, CMB constraints 
are extremely effective at constraining the available pa¬ 
rameter space. 
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